Tunneling transverse to a magnetic field, 
and how it occurs in correlated 2D electron systems 
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We investigate tunneling decay in a magnetic field. Because of broken time-reversal symmetry, the standard 
WKB technique does not apply. The decay rate and the outcoming wave packet are found from the analysis 
of the set of the particle Hamiltonian trajectories and its singularities in complex space. The results are 
applied to tunneling from a strongly correlated 2D electron system in a magnetic field parallel to the layer. 
We show in a simple model that electron correlations exponentially strongly affect the tunneling rate. 
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Tunneling plays a fundamental role in many physical 
phenomena. In the last decades much progress has been 
made in describing it beyond the one-dimensional ap- 
proximation and understanding how it occurs in many- 
body systems For charged particles, the tunneling 
rate can be conveniently controlled by a magnetic field 
applied transverse to the tunneling direction. Recently 
this effect was used to probe two-dimensional electron 
systems (2DES) in semiconductor heterostructures 
and on helium surfaces H. However, despite its interest 
and generality even the problem of single-electron tun- 
neling decay in a magnetic field (TDMF) remains largely 
unsolved. Existing results, although often highly non- 
trivial, are limited to the cases where the potential has 
either a special form P, pl||Tl"| ] (e.g., parabolic p0[), or a 
part of the potential or the magnetic field are in some 
sense weak |12|-|lq|. 

Of particular interest for the present paper is tunnel- 
ing transverse to the field from strongly correlated 2DES 
H,[17|]. In such systems exchange is small, and the tun- 
neling electron can be thought of as moving in the po- 
tential created by other electrons. This motion can expo- 
nentially strongly affect the tunneling rate. This can be 
qualitatively understood, because an electron which tun- 
nels a distance z transverse to the field B and is free to 
move in the Bxz direction, has its velocity v H in this di- 
rection incremented by lo c z (u) c — eB/m is the cyclotron 
frequency). If initially v H = 0, the energy of motion in 
the z direction is reduced by mu^z 2 /2, i.e. there arises 
a parabolic magnetic barrier for tunneling. On the other 
hand, if an electron can give the momentum in the Bxz 
direction to in-plane excitations in the electron system, 
v H can remain small, effectively reducing this barrier. 
This makes it possible to use tunneling in a magnetic 
field as a sensitive probe of in-plane electron dynamics in 
correlated systems. 

In this paper we will use an Einstein model in which the 
in-plane electron motion is a harmonic vibration about 
an equilibrium position, with one frequency (see Fig. 2 
below) [Hi . The problem is then effectively reduced to a 



single-particle problem, which mimics the many-electron 
one. As we show, the resulting tunneling exponent de- 
pends on the dimensionless parameters uj c tq and ujoTq, 
where to is the imaginary time of underbarrier motion 
for B = 0, and ojq is a characteristic frequency of in- 
plane electron vibrations. 

For smooth potentials and magnetic fields, the tun- 
neling rate can be found in the WKB approximation, in 
which the wave function is 



yj(r) = D(r) exp[iS'(r)] (h = 1). 



(1) 



Here, S(r) is the classical action. It is calculated using 
the classical equations of motion 



S = p • r, r = dH/dp, p = —dH/dr, 



(2) 



where H = (p + eA) 2 /2m + U(r) is the electron Hamil- 
tonian, and A is the vector potential of the magnetic 
field. 

In the standard approach to tunneling decay, which 
applies for B = @,|,|l9|,|o| , the action S is purely 
imaginary under the barrier. It is calculated by chang- 
ing to imaginary time and momentum in Eqs. (0), which 
then take the form of equations of classical motion in 
an inverted potential —U(r), with energy — E > —U(r). 
In the presence of a magnetic field, because of broken 
time-reversal symmetry, the replacement t — > —it, p — > 
ip, r — > r, U(r) — > — U(r) would lead to a complex 
Hamiltonian, which makes no sense and indicates that 
a more general approach is required. 

We will find the action S by solving the Hamiltonian 
equations (||) in complex time and phase space. In con- 
trast to the B = 0-case, in the presence of a magnetic 
field the action is complex for real r, i.e. the decay of the 
wave function (|l|) under the barrier is accompanied by 
spatial oscillations. The tunneling rate is determined by 
Im S at the point where the particle emerges from the 
barrier as a semiclassical wave packet, with real coordi- 
nate and real momentum. However, again in contrast to 
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the standard (B — 0) analysis, at this point the particle 
velocity is finite. Re r ^ 0. 

The trajectories (|2j) of interest for tunneling decay start 
for t — from the vicinity of the localized metastable 
state. The initial conditions can be obtained from the 
known form of the wave function ip(r) close to the po- 
tential well, both in the case where U(r) is parabolic 
near the minimum and ip is Gaussian and where 

U(r) is nonanalytic in one variable (z), which is of in- 
terest for 2DES. In both cases the trajectories (|^) are 
parametrized by two complex parameters 21,2(0), which 
can be the initial values of the in-plane coordinates 
x(0) — xi(0), y(0) = £2(0) for given z(0), and which 
in turn determine p(0) and 5(0), cf. Eq. (Q). 

To find the tunneling exponent we note that, once 
the particle has escaped, it is described by a wave beam 
which propagates in real time along a real classical trajec- 
tory r c i(t). This trajectory can be obtained by analyzing 
the fan of complex trajectories r(i), p(i) (||) for different 
£1,2(0) and finding such 21,2(0) that, for some t, both 
r(i) and p(t) become real, 



Im r(i) = Im p(t) = 0. 



(3) 



This is a set of equations for complex 21,2(0) and Im t. 
The number of equations is equal to the number of vari- 
ables, with account taken of H being real. The Re t 
remains undetermined: a change in Re t in (^|) results 
just in a shift of the particle along the classical trajec- 
tory r c i(t), see Fig. 1. Such a shift does not change Im 5. 
We note that, in contrast to what happens for B = 0, the 
classical trajectory does not have to touch the boundary 
of the classically forbidden region. 

The tunneling exponent 1Z is given by the value of Im 5 
at any point on the trajectory r c i, 



K = 2 Im 5(r c 



(4) 
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FIG. 1. (a) Complex t plane for integrating the Hamilto- 
nian equations (^) in the escape problem. The line Im t = 
const corresponds to the classical trajectory of outgoing elec- 
tron, (b) The classical trajectory r c i on the (x,z) plane. The 
bold solid lines in (a) and (b) show the "visible" part of the 
trajectory where the amplitude of the outgoing wave packet 
exceeds the amplitude of the tail of the underbarrier wave 
function. The thin solid line in (b) shows where the am- 
plitudes of the two wave functions coincide. It starts from 
the caustic (x). The data refer to the potential (^), with 
o>oTo = 1.2 and lo c to = 1.2, time in (a) is in the units of to. 



For a physically meaningful solution, Im 5 has a 
parabolic minimum at r c j as a function of the coordi- 
nates transverse to the trajectory, and the outgoing beam 
is Gaussian near the maximum. 

From (||), the tunneling exponent can be obtained by 
solving the equations of motion (Q) in imaginary time, 
with complex r. However, such solution does not give 
the wave function for real r between the well and the 
classical trajectory r c \. Neither does it tell us where the 
particle shows up on the classical trajectory. To obtain a 
complete solution of the tunneling problem, one should 
take into account the fact that 5 is a multivalued func- 
tion of r, even though it is a single-valued function of t 
and 21,2(0), i.e. several trajectories (||) with different t 
and 21,2(0) can go through one and the same point r. 
The wave function is determined generally by one of the 
branches of 5(r). 

Branching and multivaluedness of the action are famil- 
iar from ID tunneling problem, where S — St cx (z — z t ) 3 ^ 2 
near the turning point z t p9[. In multidimensional sys- 
tems, branching generally occurs on caustics I n our 
problem, in contrast to the usually considered case, the 
trajectories r(t) will be complex, as will also be the caus- 
tics. Caustics of most general type are envelopes of the 
trajectories r(t) (0). They are given by the equation 



J(r) = 0, J(r) 



d{x,y,z) 



9(2i(0), x 2 (0),t) 



(5) 



The prefactor in the WKB wave function (|l]) is D = 
const x J~ x / 2 . Therefore the WKB approximation does 
not apply close to the caustic (cf. [pl| ). 

The caustic of interest is the one where the analyti- 
cally continued wave functions of the outgoing semiclas- 
sical wave and the WKB tail of the intrawell state are 
connected. The amplitude of the semiclassical wave in- 
cident on the barrier from z — > 00 should be set equal to 
zero. Local analysis near the caustic is similar to that in 
the ID case. It is convenient to change to the variables 
2', y' , and z' which are locally parallel and perpendicular 
to the caustic surface, respectively. We set z' = on the 
caustic. For small \z'\, 



5(2', y', z') w 5(2', y', 0) + aiz' + a 2 z' z ' 2 



(6) 



(the coefficients Q4.2 = o-i^x' ,y') can be expressed in 
terms of the derivatives of 5, r over 21,2(0), t on the caus- 
tic). As in ID problem, the prefactor in p) depends on 
the distance to the caustic as Doc (z') -1 / . However, in 
the present case 5 contains a linear term aiz', and there- 
fore the momentum perpendicular to the caustic is finite. 
We note that Eq. (||) and the condition Im r = define 
a line in real space, which can be called a caustic line. 

For real r, there is a switching surface which separates 
the ranges where Im 5(r) is smaller for one or the other 
of the solutions connected on the caustic [^2j . Only the 
solution with the smaller Im 5 should be held in the 
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WKB approximation. It is this condition that determines 
where the outgoing wave shows up from beneath the tail 
of the intrawell state. The switching surface starts from 
the caustic line, where the branches of S merge together. 
The cross sections of the switching surface and the caustic 
line by the plane (x, z) are shown in Fig. 1. 

We now apply these general results to a simple model 
which is relevant to electrons on helium. The correspond- 
ing geometry is shown in Fig. 2a. For typical densities 
n ~ 10 s cm~ 2 and temperatures T < 1 K, these electrons 
form a Wigner crystal [^3| or a nondegenerate liquid Q . 
In both these cases, the characteristic frequency of vi- 
brations about a (quasi) equilibrium in-plane position is 

co = (2^e 2 n 3 / 2 /™) 1/2 . 

The initial conditions for the equations of underbarrier 
motion (|^) can be chosen at an arbitrary plane z = const 
close to the electron layer and yet deep enough under the 
barrier so that the electron wave function is semiclassical. 
We set z = z(Q) = on this plane. In the Einstein 
model, the potential for in-plane motion for z = is 
(muj 2 /2)(x 2 + y 2 ), and ij){x, y,z — 0) is Gaussian in x, y. 
Then for t = in (§), 

z(0) = 0, Pz (0) = i 7l S(0) = [x 2 (0) + y 2 (0)] , (7) 

where 7 = [2mU(r — 0)] 1 / 2 is the reciprocal localiza- 
tion length in the z-direction, 7 3> (mujo) 1 ^ 2 (the en- 
ergy of the localized state is set equal to zero). The 
initial values of the in-plane momentum components 
Pj(0) = dS(0)/dxj(0) — imoJoXj(0) (where j = 1,2 enu- 
merates the in-plane coordinates). 

For electrons on helium, tunneling occurs if there is 
applied an electric field £±_ which pulls electrons away 
from the helium surface j2j|. The tunneling barrier is 
formed by the image potential and the potential of this 
field. It is nearly triangular at distances from the surface 
which are much greater than the effective Bohr radius 
7~ 1 . The barrier width L = j 2 /2me£± for B = 0. Typ- 
ically L <C n -1 / 2 , and therefore the effective electron 
potential energy is well represented by 

a(r) = =V + „»> + i(l-±) (*>»>■ («) 

With (j^), the equations of motion (g) become linear 
and can be readily solved. The symmetry of the po- 
tential U(x,y,z) = U(±x,±y,z) gives rise to a specific 
symmetry of the set of the trajectories (||), and one can 
show that the caustic of interest intersects the real space 
for x = and some z = z c (z c = L for B = 0). 

For B along the y-axis (see Fig. 2a), the motion in 
the y direction is decoupled and the problem becomes 
two-dimensional. For z < z c , the function Im S has two 
branches each of which is symmetrical in x. The branch 
1 describes the tail of the intrawell wave function before 
branching. It has a minimum at x = for given z, and 



monotonically increases with \x\ and z. As expected, the 
slope d Im S/dz is finite at the branching point z c . The 
branch 2 corresponds to the wave "reflected" from the 
caustic. This branch is nonmonotonic in z for x = 0, 




FIG. 2. Two branches of the action on the symmetry axis 
x = as a function of the tunneling coordinate z before the 
branching point, for the same parameters as in Fig. I. The 
vicinity of the cusp at z c is zoomed in inset (b). Inset (a): 
the geometry of tunneling from a correlated 2DES transverse 
to a magnetic field; ujq is the Einstein vibration frequency in 
the Wigner solid. 

with a minimum at z m < z c . For z m < z < z c , it has 
two symmetrical minima for x ^ 0. They lie on the 
classical trajectory shown in Fig. lb, and merge together 
for z = z m . 

As discussed above, Im S is constant on the classi- 
cal trajectory in Fig. lb. This trajectory goes through 
the point x = 0, z = z m and is symmetrical in x. Al- 
though the potential U (r) is even in x and is minimal for 
x = 0, the escaped particle "shows up" on the classical 
trajectory for finite x. This happens where Im Si(r c \) = 
Im 52(r c i) =72-/2 (the subscript enumerates the branches 
in Fig. 2). The particle has finite velocity and moves away 
from the barrier. 

Since the point x — 0, z — z m lies on the classical 
trajectory of interest (although on the section "hidden" 
by the tail of the intrawell state), the tunneling exponent 
is given by 1Z — 2ImS , 2(x = y = 0,z m ) and can be 
calculated in imaginary time, with imaginary x(0) 

K= -is 2 t 3 - 3vo(l - t) 2 + 3(1 + m)(l - r) (9) 
+iv 2 r, K = 2-fLK/3v 2 , 

where vq = u)qTq, v c = lo c tq are the dimensionless in- 
plane and cyclotron frequencies (tq — 2mL/"f is the 
"duration" of underbarrier motion in imaginary time for 
B = 0), v 2 = v 2 + v 2 , and r = it /to is given by the 
equation 

[v 2 vq{\ - t) - v 2 \ tanh^r = v{v%t - v 2 ). (10) 

The tunneling exponent as a function of luq, lu c is shown 
in Fig. 3. For loq — (no electron-electron interac- 
tion), the magnetic barrier makes tunneling impossible 
for uj c tq >1 1 (see curve 1; r — ► 00 for vq — Q^v c ^ 1). 
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Even comparatively weak in-plane confinement elimi- 
nates this effect. The reduction of the tunneling suppres- 
sion is significant already for small uiqTq, and increases 
fast with increasing wo T o- 
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FIG. 3. The dependence of the tunneling rate on mag- 
netic field, W = W(B)/W(0). The curves 1 to 4 re- 
fer to ojoto = 0, 0.2, 0.4, 0.6. Magnetic field eliminates 
single-electron tunneling for u> c tq > 1 (cf. curve 1). 
Inset: tunneling exponent vs in-plane frequency ljq f° r 
lOcTq = 1.0,2.0,3.0 (curves a,b,c). 

The above results provide an explanation of the mag- 
netic field dependence of the tunneling exponent for elec- 
trons on helium, which was observed to be much weaker 
than it would be expected from the single-electron 
theory. Detailed comparison with the data H will be 
discussed elsewhere where the model will also be 
extended in order to include the realistic vibrational spec- 
trum of the Wigner solid. At zero temperature this ex- 
tension does not change the results significantly, because 
electron tunneling is accompanied by excitation of mostly 
short-wavelength vibrations, which are reasonably well 
described by the Einstein model used above. 

In conclusion, we have shown that, under suitable con- 
ditions (lo c to <^ 1, luoTq > 1), correlations in a 2DES can 
exponentially strongly affect the rate of tunneling escape 
transverse to a magnetic field. We have also shown that 
the problem of single particle tunneling in a magnetic 
field can be solved in the semiclassical limit by analyzing 
the Hamiltonian trajectories of the particle in complex 
space and time. The connection of decaying and prop- 
agating waves occurs on caustics of the set of these tra- 
jectories. This approach does not require us to consider 
any piece of the electron potential or the magnetic field 
as a perturbation. It gives us an escape rate which is 
generally exponentially smaller than the probability for a 
particle to reach the boundary of the classically accessi- 
ble range U(r) — E. Finally we have obtained explicit 
results for a simple model of an electron tunneling from 
a helium surface transverse to a magnetic field. 
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